Global selection rule in chemical coupling 



< 



X 



2 



Asoke P. Chattopadhyay 1 and Anjan K. Dasgupta 
1 Department of Chemistry, University of Kalyani 

bJQ' 

s : 2 



Kalyani 741235, India 
Department of Biochemistry, University of Calcutta 



O ' 35 Ballygunje Circular Road, Calcutta 700019, India 

m 

S: 

% ■ Abstract 
q 

■ Coupling and decoupling of chemical reactions are explored 

•l-H . 

through a modified heat balance equation. Reaction enthalpies 
are found to play crucial role; the sign of their product for a pair 
of consecutive chemical reactions determine whether they couple 
£^ ■ or not. The possibility of a coupling-uncoupling transition for 

i> . 

such reactions is thus introduced for the first time. The present 

OO 

work resolves a paradox concerning negative efficiency of coupled 

■ consecutive chemical reactions. Enthalpy is also shown to be a 
*5a ■ "constant of motion" along the reaction coordinate as long as the 

■ mass action ratio varies little with temperature. The present anal- 
ysis puts the observed difference between calorimetric and van't 
Hoff enthalpies on a quantitative basis. A case study is presented 
with a third order reaction where the enthalpic criterion is useful 
in selecting the more probable of the alternative mechanisms. 
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I. INTRODUCTION 



The present work extends an earlier attempt M to generalize the conven- 
tional phenomenology for describing thermodynamics of chemical reac- 
tions. The conventional approach, though based on the limiting assumption of 
near equilibrium, was applied to complex processes 0] successfully. Coupling 
of chemical reactions, however was treated [Q using a rather special example 
of a cyclic reaction system. For non-cyclic e.g. consecutive reaction schemes, 
the Onsager matrix is shown to be diagonal, ruling out phenomenological cou- 
pling. Conventional chemical wisdom, on the other hand, assumes a set of 
reactions to be coupled provided they have common intermediate(s). More 
importantly, the efficiency of coupling can be shown to be always negative for 
coupling schemes treated in the traditional way, whether in linear or n non- 
linear domain. This renders the description of coupling itself as invalid in the 
established format. We show, on the contrary, how minor modifications in the 
balance equations lead to alteration in the definition of macroscopic fluxes and 
forces in chemical reaction systems |||7]] permitting coupling. The reactions 
may now be coupled if enthalpy change in each is non-vanishing. Also, en- 
thalpy remains a " constant of motion" along reaction coordinate provided the 
mass action ratio has a small variation with temperature. The recent debate 
on difference between van't Hoff and calorimetric enthalpies || are put on a 
more quantitative basis with our approach. A case study with a third order 
reaction having two possible mechanisms shows that our approach correctly 
identifies the more probable pathway. 
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II. THE BALANCE EQUATION REVISITED 



Let us briefly recall the phenomenological equations for chemical coupling 
widely used in chemical literature. j|,[7| For a set of reactions p(— 1, ....r) in a 
fixed volume, the rates of change of entropy, concentration etc. are given by 

ds v /dt = -V.J s + cr (1) 

dCi/dt=-V.Ji + ^ipVp (2) 
p 

dq v /dt = -V.J, + £ v p AH p (3) 

Tds„ = dq v - £ /^dc* (4) 

i 

are the density and flux terms for entropy, Cj, Jj those for the ith 
chemical species and q v , J q those for heat. V{ p is the stoichiometric coefficient 
of the i th species in the p th chemical reaction for which the velocity is v p . a 
is the entropy production term. While AH P is the enthalpy change of the p th 
chemical reaction, the corresponding Gibbs free energy change, AG P , is given 
by the law of mass action as 

AG P = -RTln(v+/v-) (5) 

where if}" and v~ are the forward and reverse reaction velocities of the p th 
recation, and v p = v+ — v~ . From eqs. (1) - (4), equating non-gradient parts, 

a = J q V(l/T) - Ji-V(JH/T) + J2 v p( AH p ~ Y,Wi P )l T ( 6 ) 

i Pi 
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For an isothermal chemical reaction system in a well-stirred (or homoge- 
neous) medium we get, 



a = 



Y, v p as p 



(7) 



p 



from AG p = AH p — TAS P and the second law. Note that the rate of entropy 
production is obtained as a stoichiometric sum of entropy changes of reaction 
steps, without invoking any assumption of linearity of processes. 

Our eq. (7) can be compared with the standard one for a found in chemical 
literature, j|,[7j viz. 



To 




p 



Following standard phenomenological notations, 



£^pp' AG p 



(9) 



p 



Linearising eq. (5), and using v p = v p — v with eq. (9), we obtain 



L / 

pp 



(10) 



where 5 pp / is the Kronecker delta. Such a diagonal nature of L makes coupling 
between two different reactions impossible. Please note that coupling between 
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cyclic reactions can still emerge in this treatment, Onsager's example [|J being 
the most famous one. 

This impossibility does not occur in our treatment since unlike eq. (8), eq. 
(7) leads to a different phenomeno logical equation for the reaction velocity 

p 

III. PHENOMENOLOGICAL VS. CHEMICAL COUPLING 

Coupling of chemical reactions, while gaining in popularity over the past 
few decades, have received little attention from theoreticians. The major con- 
tributors have been Prigogine and his coworkers |4],[J . Perhaps, the importance 
of coupled reactions is felt nowhere more than in treatments of biochemical 
cycles ||10|| . Glycolytic or the basic bioenergetic cycle (oxidative phosphoryla- 
tion) are examples of intricate coupling of consecutive and cyclic reactions |J. 
Yet, quantitative expression of reaction coupling is absent in existing literature 
in these fields. 

According to Prigogine @], a pair of reactions with affinities A\ and Aj and 
velocities V{ and Vj can couple if signs of A^Vj and AjVj are different. This 
means that one reaction must drive another for the two reactions to couple. 
Based on this idea, Rottenberg || defined efficiency of coupling as 

7] = -A 1 v 1 /A 2 v 2 

where the subscripts 1 and 2 denote the driven and the driving reactions re- 
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spectively. However, Hill ]TT|] pointed out that 



A x v > 

for individual reactions. Thus efficiency of coupling is negative and coupling of 
consecutive reactions at least is ruled out. This poses a far stronger challenge 
to the conventional theory describing coupling, as an objection to the problem 
posed by diagonal nature of the Onsager matrix, as shown in eq. (10), could 
be that the treatment there is strictly linear. Although Hill's derivation is 
consistent with eq. (5), his conjecture that coupling of chemical reactions is 
only possible via common intermediates and has no phenomenological meaning 
otherwise, is hardly acceptable. In fact, Hill never sought to explain why a 
consecutive reaction pair always has a negative efficiency. The lacuna was in 
not appreciating a more fundamental problem, viz. the nature of the driving 
force behind reaction coupling. This impasse has been resolved by the present 
work, where the key steps in this regard were modification of eq. (8) to (7), 
and of eq. (9) to (11). From eq. (11) it is clear that in our theory the efficiency 
of coupling need not be always negative. Also, the Onsager matrix L / is not 
necessarily diagonal. It can thus be used as a proper measure of the coupling 
strength between two reactions. The general approach is not only useful in 
resolving the paradox stated earlier, it also provides a simple thermodynamic 
criterion for coupling in simple reaction systems. For example, if the enthalpy 
of any reaction step p is zero, that reaction can not couple with any other 
reaction step p as L pp > again becomes diagonal. 
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IV. GIBBS RELATION ALONG REACTION TRAJECTORY: AH C al 

VS AH VH 

This is an interface between kinetics and thermodynamics of reaction pro- 
cesses 0. For a reaction system the free energy change of the p th reaction 
away from equilibrium can be expressed by 

AG P = AG° p + RTlnKp (12) 
where the mass action ratio K p is given by 

K P = E C ? P (13) 

3 

Writing the entropy change along the p th reaction as a stoichiometric sum of 
molal entropies, Si 

ASp = E*v^ (14) 

i 

From the relation AS P = -dAG p /dT, eq. (12) and AG P = AH P - TAS P , 

AS P = AS° p - RlnK p - R(dlnK p /dlnT) (15) 
A simple consequence of eq. (15) is, 
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AH P = AH° p - RT{dlnK p /dlnT) 



(16) 



Let us identify the two AH terms appearing in eq. (16) clearly. At equilib- 
rium, eq. (16) becomes 

AH e p q = AH° p - RT[dlnK p /dlnT} eq (17) 

If the second term on RHS becomes very small, we are left with 

AH P ~ AH° p (18) 

Eq. (18) may be stated as follows: the enthalpy remains approximately a 
constant of motion along the reaction co-ordinate for any reaction (the van't 
Hoff approximation). Note that where it is not so, i.e. where AH p depends 
on the reaction coordinate and may even change its sign from AH p , a pair of 
coupled reactions may become decoupled or vice versa, as the coupling depends 
on sign of the product of the enthalpies of the respective reactions. A coupling- 
uncoupling transition is therefore possible for such a set of chemical reactions. 
Much further work is needed in this area. 

AH p q is clearly the experimentally measured enthalpy change in a reaction 
viz. AH cai, which is conventionally measured at equilibrium (or near equi- 
librium) conditions. But the first term on RHS of eq. (16) or (17) can be 
identified with the van't Hoff enthalpy of the reaction viz. AH vH , defined as 
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AH vH = AH° p = RT(dlnK eq /dlnT) 



(19) 



Based on this discussion and the last two equations, we can write 



AH cal = AH vH - RT[dlnK p /dlnT] eq (20) 



There is a large and growing body of evidence |8|]i2Hl9l (including some un- 



published work ||20[| ) of discrepancy between AH ca i and AH v h- There is also 



an awareness of the importance of the ratio of these two values, especially in 



interpreting biocalorimetric data ||21|| . While for simple chemical reactions the 
ratio is close to unity, for reactions involving macromolecules e.g. in protein 
folding, there is clear departure of the ratio from unity. In biochemical lit- 



erature, |19| , |2~TH the numerical value of the ratio (which may vary from 0.5 to 
more than 4 ||13|| , say) is taken to provide a measure of cooperativity of the 
biochemical reaction (e.g. folding). Our analysis gives a clear insight into 
the difference between the two enthalpy values. The origin of this difference 
stems from the difference in temperature dependence of equilibrium and non- 
equilibrium mass action ratios. Let us also recall that such difference between 
equilibrium and non-equilibrium values are accepted naturally for Gibbs free 
energy, for example, and the difference considered in terms of a mass action 
ratio. 

We realize that further simplification of eq. (20) may be difficult. Instead, 
we use certain experimental data to show that the mass action ratio, K p , may 
have a scaling dependence on the experimental temperature. Holtzer ]12| es- 
timates that the difference between the calorimetric and van't Hoff enthalpies 
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may be of the order of 45 calories per stoichiometric unit for simple chemical 
reactions. From eq. (20), this leads to 

dlnK p /dlnT ~ -45/(300 x 1.98) = -0.075 (21) 

at room temperature (300° K). We immediately obtain 

K p ~ const x T~ a (22) 
where a ~ 0.075 for simple chemical reactions. 

V. COUPLING COEFFICIENTS IN TWO STEP REACTION - 
KINETIC APPROXIMATION 

Using the dissipation eq. (7), we may express the phenomenological rela- 
tion of a pair of reactions by 



AS 2 = R21V1 + R22V2 



where V\ and t>2 are velocities of the two reactions and R^ are the coupling 



terms. The pre-equilibrium condition [22] is given by 



Ad 



and V\ 



(23) 



If we use the approximate linear relation in such cases 
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AS! = R U V! + #12^2 ^ #12^2 (24) 

AG i ~ implies AHi ~ TAS±. Using the arguments of eq. (24), the entropy 
change in the first reaction is given as 

AS i ~ AHl/T (25) 

Using eq. (24), eq. (25) assumes the form 

AH® = R u v 2 (26) 

where R l2 = TR 12 . The approximate form of the dissipation equation can be 
expressed as 

a = ASiUi + AS 2 v 2 ~ A52V2 > (27) 

As v 2 , velocity of the rate-determining step, is positive both for positive and 
negative coupling, the pre-equilibrium approximation should satisfy 

AS 2 > (28) 

Eqs. (25) - (28) express the nature of thermodynamic coupling that may 
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exist for kinetic schemes satisfying the pre-equilibrium condition. Eq. (26) 
implies that velocity of the rate determining step will be proportional to the 
enthalpy change in the pre-equilibrium step. It may be interesting to note that 
for positive coupling, the pre-equilibrium step must be endothermic and vice 
versa. If on the other hand the coupling is negative the relation AH® < must 
be satisfied. The negative coupling implies that one of the reactions among 
the pair has a negative dissipative component and is therefore driven by the 
other. 



VI. A CASE STUDY WITH IN A SIMPLE REACTION 



Let us consider a typical chemical reaction 



2NO + 2 ^ 2N0 2 



Using reported molar enthalpy values |23] the overall AH of this reaction 
is approximately 27.02 Kcal/mole. The standard mechanism (henceforth re- 
ferred to as mechanism I) cited p4| is the following one: 



2NO = N 2 2 
N 2 2 + 2 ^ 2N0 2 



This mechanism is able to explain the accepted rate law for the overall reaction 
viz. 
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Rate = k[NO} 2 [0 2 } 



as also the negative sign of the activation energy of the overall reaction. [^5| 
Recently, |26| an alternative mechanism (referred to as mechanism II), has 
been proposed, namely 



NO + 2 = OONO 
OONO + NO^ 2N0 2 



which agrees with the rate law given above. According to Plambeck, [26 spec- 
troscopic evidence admits of simultaneous existence of both these mechanisms 
in the overall reaction, although concentration of OONO may be larger than 
N 2 2 . Both the intermediates are short lived species, and although there is 



speculation about the structure of N 2 2 , [27 nothing is known about OONO. 

We performed ab-initio calculations on N 2 2 and OONO. The GAMESS 
software package was used with the ccpVTZ basis set, [^J augmented by a 
d and diffuse s, p type polarization functions. We present only thermochemical 
data in Kcal/mol (for AH or AG ) or cal/mol-K (for AS" ). Subscripts 1 and 
2 represent the two steps in the overall reaction. Along with standard data 
|23H27| our results are (with an estimated error limit of ±10 percent) 
For mechanism I: 



AS i = -40.68 cal/mol-K 
AH° = -30.2 Kcal/mol 
AG? = -17.9 Kcal/mol 



AS 2 ° = 5.93 cal/mol-K 
AH$ = 3.18 Kcal/mol 
AG° 2 = 1.4 Kcal/mol 
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Corresponding values for Mechanism II are 



AS? = -34.34 cal/mol-K 
AH° = -11.60 Kcal/mol 
AG° = -1.29 Kcal/mol 



AS 2 ° = -0.39 cal/mol-k 
AH$ = -15 Kcal/mol 
AG° 2 = -15.29 Kcal/mol 



It is clear that in both the mechanisms the enthalpy change in the pre- 
equilibrium step has a negative sign. Eq. (26) implies that one reaction 
must drive the other. From eq. (28) it follows that for both the mechanisms, 
AS 2 = AS® — RlnK 2 > 0, where K 2 is the mass action ratio for the second 
step of the reaction. Therefore, 
for Mechanism I 



If the intermediate concentration is small, AS 2 must be greater than a pos- 
itive quantity. This need not be the case if the intermediate has a higher 
concentration, since then the right hand side of inequality (30) shifts towards 
more negative value. Incidentally, the first mechanism has a positive AS 2 . In 
the second mechanism AS 2 has a value approaching zero. In this case, (30) 



AS° 2 > Rln[N0 2 } 2 /([N 2 2 }[0 2 }) 



(29) 



and for Mechanism II 



AS 2 ° > Rln[N0 2 } 2 /{[NOf[OONO}) 



(30) 
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can hold good provided the intermediate concentration has a higher value. As 
pointed out by Plambeck, |2C| both the mechanisms are known to exist, but 
higher concentration is seen for OONO. The observation is thus in accordance 
with the thermodynamic arguments presented above. 

VII. CONCLUDING REMARKS 

The present paper shows that coupling of chemical reactions should not 
be described in vague qualitative terms e.g. existence of common chemical 
intermediates. A pair of reactions remain decoupled as long as the product of 
their reaction enthalpies is non-negative. Coupling can only occur if they have 
enthalpies of opposite signs. An interesting corollary is that the same pair of 
reactions can undergo a transition from coupled to decoupled state (or vice 
versa), provided the mass action ratio of one or both the reactions change ap- 
preciably with progress of the reaction(s). An important modification brought 
forward by the present work is that the rate of internal entropy change near 
equilibrium is a weighted sum of the entropies (and not free energies or chem- 
ical affinities, as the form in which it is usually expressed) of the participating 
reactions. The observed deviation of calorimetric (measured) enthalpy from 
the van't Hoff value for a reaction has also been explained, the origin of which 
is shown to be in the difference in the temperature dependence of equilibrium 
and non-equilibrium mass action ratios. For a class of simple chemical reac- 
tions, where the deviation is small, the enthalpy remains a constant of motion 
along the reaction coordinate. Finally, our treatment is shown to identify the 
more probable of alternate pathways for a typical third order chemical reac- 
tion. 
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